Nonlinear a model for disordered superconductors 



I. V. Yurkevich, and Igor V. Lerner 
School of Physics and Astronomy, University of Birmingham, Edgbaston, Birmingham B15 2TT, UK 

(24 January 2001) 

We suggest new variant of the nonlinear a-model for the description of disordered superconduc- 
tors. The main distinction from existing models lies in the fact that the saddle point equation is 
solved non-perturbatively in the superconducting pairing field. It allows one to use the model both 
in the vicinity of the metal-superconductor transition and well below its critical temperature with 
full account for the self-consistency conditions. We show that the model reproduces a set of known 
results in different limiting cases, and apply it for a self-consistent description of the proximity effect 
at the superconductor-metal interface. 



I. INTRODUCTION 

Since a seminal paper by Wegner,0 a field theoretic 
approach to disordered systems based on the nonlinear 
a model (NLcrM) became one of the most powerful tools 
in describing localisation effects and mesoscopic fluctua- 
tions. The main advantage of this approach lies in for- 
mulating the theory in terms of low lying excitations 
(diffusion modes) which greatly simplifies perturbative 
and renormalization group calculations, and on the other 
hand allows a non-perturbative treatment. 

Such an approach has been successfully extended 
to the description of disordered superconduciors.Bo It 
was based on the Fermionic representationO of Weg- 
ner's NLcrM extended to include the electron-electron 
interactional. The starting point in these workso □ was 
a microscopic model of interacting electrons in a ran- 
dom potential. The effective NLcrM includes an extra 
Bosonic field describing the superconducting order pa- 
rameter A. Then the lowest-order expansion in A is 
used. This makes such an approach a good working tool 
in the vicinity of the superconducting transition where 
all the interaction channels can be easily included which 
makes it very useful in describing different aspects of the 
metal-superconducting transitions. 

An alternative r-approach to the NLcrM for dirty 
superconductorsQ~tL2l starts from the Bogoliubov-de 
Gennes equations (or, equivalently, Gorkov's equations) 
without imposing a self-consistency condition on the 
superconducting order parameter A which is consid- 
ered as given. Then the initial many-body problem 
turns into a single particle one which makes applica- 
ble powerful techniques based on the supersymmetric 
NLcrM .EJ Such a supersymmetric approach has been re- 
cently developed in Refs. |l^ and applied to the descrip- 
tion of non-perturbative aspects of the proximity effect 
in superconducting-normal-metal structures. In this ap- 
proach A was taken into account just by the boundary 
conditions (Andreev reflection) for the normal region. A 
natural disadvantage of this (and any supersymmetric) 
approach is that no interaction can be included beyond 



the mean-field approximation; thus it is impossible to 
describe an effect on the superconducting order param- 
eter of disorder in the normal metal (or even inside the 
superconducting region) . 

A novel NLcrM developed in this paper starts from 
a microscopic model of electrons in a random potential 
with the BCS attraction, and the order parameter A is 
treated as a dynamical field, similar to the earlier devel- 
oped microscopic approachxTu We are using the standard 
fermionic replica approachO in temperature techniques 
For a long time, it was widely believed that such an 
approach cannot be used for nom-ipprturbative analysis. 
However, it was recently shownti'O that this is not the 
case, since the well-known exact non-perturbative result 
was reproduced fxom the fermionic replica NLcrM, as well 
as more recentlytj within the Keldysh technique. 

In the initial approachQ to interactions within the 
NLcrM, a saddle point approximation was identical to 
that of the non-interacting— problem. This scheme was 
recently greatly improvedE3 by choosing (within the 
Keldysh technique) the saddle point with taking account 
of the interaction which considerably simplified any fur- 
ther analysis. Such an analysis has been directly ex- 
tended to dirty superconductors in Ref. We consider 
a model where, for simplicity, the Coulomb repulsion is 
not included. A distinctive feature of our approach is 
a change of the saddle point (and of a subsequent ini- 
tial approximation) in the presence of the superconduct- 
ing order parameter. This is similar but not identical 
to the choice suggested in Ref. [l^ (when applied to the 
Coulomb interaction, it would lead to a different variant 
of the NLcrM). The NLcrMlij is optimized to maximally 
simplify the lowest perturbational order while by sacri- 
ficing this we arrive at quite a general formulation of the 
model with different specific approximations being made 
for different applications. 

As usual, we restrict our consideration to the limit 
of dirty superconductors when A <C l/r c i <C ef (or, 
equivalently, vpT c \ <C £ where r i is the elastic mean free 
time, and £ in the correlation length in dirty supercon- 
ductors). After describing describing in detail an alterna- 
tive saddle-point approximation, we show how the model 
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reproduces a set of known results in different limiting 
cases, and apply it for a self-consistent description of the 
proximity effect at the superconductor-metal interface. 



II. BASIC MODEL 

We consider the standard BCS Hamiltonian in the 
presence of a random potential, u(r). For complete- 
ness, we start with outlining the standard procedureo of 
a field theoretic representation in the temperature tech- 
nique for this Hamiltonian. The corresponding action 
has the form 



S — Sq + Si 
S = dxtp*(x) 







— + £ + u(r) 

OT 
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Here ip a (x) is a Grassmannian fieloQi anti-periodic in 
imaginary time r with period 1/T, x = (r,r), s = (|, |) 
is the spin index, Ao is the BCS coupling constant, and 
from now on we set h = 1. 

The random potential u(r) is supposed to be Gaussian 
with zero mean and the standard pair correlator, 



<«(r)u(r')> = 



1 



I'KVT, 



cl 



*(r-r'), 



(2) 



with v being density of states, and r e \ elastic mean free 
time. The operator £ in (fl^) is defined as 



2m V 



■tV 
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where A is a vector potential of an external magnetic 
field. 

Averaging over u with the help of the standard replica 
trick gives the quartic in ip term in the action. Using 
the Hubbard-Stratonovich transformation, one decouples 
both this term and the BCS term, Eq. ([k]), the former 
with the help of a matrix field a = a(r;r,T') and the 
latter with the help of a pairing field A = A(r; r) which 
will eventually play the role of the order parameter. This 
results in the following effective action: 



1 J 8r c i 



+ / dx^(x) 
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dx|A(a;)| 2 + 
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Here the replicated Grassmannian fields are 



*^(C*) T = ^=(^,-^), 



where i — 1 . . . N are the replica indices (N = in 
the final results). The standard doubling of these fields 



(ip — > ty) is convenient to separate diffuson and cooperon 
channels for electrons propagating in the random poten- 
tial; C is the charge conjugating matrix defined by the 
above equation. The matrix fields a and A are defined 
in the space spanned by \& ® ^ which is convenient to 
think of as a direct product of the N x TV replica sector, 
2x2 spin sector, and 2x2 'time-reversal' sector. The 
field a is defined by its symmetries, 



Ca T C- 



(4) 



and Tr in Eq. (g) refers to a summation over all the ma- 
trix indices, an integration over r and a double integra- 
tion over r (as a is not diagonal in t). 

The field A is an Hermitian and self-charge-conjugate 
matrix field, which is diagonal in the replica indices and 
coordinates r and r, and has the following structure in 
the spin and time-reversal space: 



A = - ( A' • 



A'" 



(5) 



where A' and A" are real and imaginary parts of the 
(scalar) pairing field A; f* r and are Pauli matrices 
(i = 0,1,2,3 with f = 1) that span the time-reversal 
and spin sectors, respectively. 

The integral over electron degrees of freedom is per- 
formed in a usual way, so that one reduces the effective 
action (in the Matsubara-frequencies representation) to 
the following form: 



-Tr In 



2r n 



(6) 



Here e = diage„, while e n — n{2n + 1)T is the Fermionic 
frequency, and uj = e — e' is the Bosonic one. Since A is 
diagonal in the imaginary time r, it is a matrix field in 
the Matsubara frequencies. 

The action is a standard starting point for a further 
field-theoretic analysis. To construct a working model, 
one needs to expand in some way the Tr In term in 
Eq. (||) . Our goal here is to derive a field theoretic model 
which is fully self-consistent in terms of the supercon- 
ducting order parameter A and does not use a small-A 
expansion. We restrict our considerations to the limit of 
dirty superconductors when A <C 1/t c \ <C £f- Other- 
wise, we do not impose any limitations on A, and will 
derive the model applicable both in the vicinity of the 
transition and deeply in the superconducting regime. 
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III. SADDLE POINT 



Our starting point is to construct a saddle point ap- 
proximation to the action (Q) in the presence of the field 
A. As usual, we vary the action with respect to the field 
a which gives 



(J) 



As l/r e i is much greater than both temperature T and 
the order parameter A, the matrix e + A plays the role 
of a symmetry breaking field. We look for a solution in 
a way similar to that in the metallic phase where such 
a role is played by the matrix e alone. In the metallic 
phase, the saddle-point equation with e ^ has a unique 
solution <t = A, where A is diagonal in e and unit in the 
replica and spin sectors: 



A = diagjsgne} , 



(8) 



For e = 0a degenerate solution to the saddle point equa- 
tion is given by any matrix a of the symmetry (^) obeying 
the condition a 2 = 1. Such a matrix can be represented 
as a = [/^ALL,with U belonging to an appropriate sym- 
metry group S3 

Similarly, a solution to Eq. (Q) in the presence of e + A 
is given by 



&,. p . = VlAV A , 



(9) 



where Va is the matrix that simultaneously diagonalises 
both a and e + A. This means that it should be found 
together with yet unknown eigenvalues A = diag A e from 



A 



A ■ 



(10) 



Naturally, one expects Va to become unit matrix above 
the superconducting transition temperature T c . 

Assuming that both the fields A(r) and cr(r) are 
smooth functions of r and looking for a spatially inde- 
pendent solution to Eq. (0) (i.e. ignoring at this stage 
the fact that £ and Va do not commute), one substitutes 
expressions (^|) and (jl(]) into Eq. (^) thus reducing it to 



<+— A + iX 



The scale of A is defined by e ~ T and A which are both 
■C 1/t c i in a dirty superconductor. Thus it is easy to 
verify that the saddle point is given by Eq. (^) with the 
eigenvalues A, Eq. (||), being not affected by the presence 
of superconductivity. Let us stress that this saddle point 
is obtained by a non-perturbative in A rotation (^) of the 
metallic saddle point A. This should lead to an effective 
functional valid anywhere in the superconducting phase 
rather than only in the vicinity of T c . 



Such an effective functional which includes fluctuations 
around the saddle point is obtained in the standard way. 
First, one constructs a saddle-point manifold of matri- 
ces a obeying the saddle-point equation at A = 0, and 
then one expands the Trln term in Eq. (|) in both the 
symmetry breaking term A and gradients of the fields V. 
The saddle-point manifold is convenient to represent as 
follows 



o = VlQV A 



Q = U^AU, 



(12) 



where Q represents the saddle-point manifold in the 
metallic phase and a is obtained from Q by the same 
rotation (^) as a s .p. is obtained from the metallic sad- 
dle point A. Therefore, Q is defined, as in the metallic 
phase, on the coset space S(2N)/S{N) ® S(N) where, 
depending on the symmetry, S represents the unitary, 
orthogonal or symplectic group. Before describing the 
expansion, let us stress that one could expand the Trln 
term without making the rotation (|l2|), i.e. in powers 
of Vct and of e + A. Although this would be formally 
an expansion within the same manifold, performing first 
the rotation (Jl^) simplifies enormously all the subsequent 
considerations and leads to a new variant of the nonlinear 
cr model. 

After substituting Eq. ( |l2] ) into Eq. @, one obtains 
the following representation for the Tr In term: 



SS = -iTrlnjGo- 1 + V A [i,v£\ - t(UXU^)} , 



where 



The expansion to the lowest powers of gradients and A is 
easily performed and results after some straightforward 
calculations in the following action: 



5 = ^£ 



dr|A.| 2 + ^Tr 



j (dQf - XQ 



(13) 



where Tr refers to a summation over all the matrix indices 
and Matsubara frequencies, as well as to an integration 
over r. The long derivative in Eq. (13) is defined as 



dQ = VQ 



le 
c 



A A --Ar, tr , 



Q 



d Q + [A A ,Q] , (14) 



(11) where the matrix A a is given by 



A A = V A d Q vl , 



(15) 



and do is the long derivative (|1J) in the absence of the 
pairing field A. Both V A and A should be found from 
the diagonalization of e + A, Eq. (|l0|). Although such 
a diagonalization cannot be done in general, it will be 
straightforward in many important limiting cases. For 
A = 0, the field Aa vanishes, d — > do, and A —> e, so 
that the functional ( p"3| ) goes over to that of the standard 
nonlinear a model for non-interacting electrons. 
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The cr model defined by Eqs. ©-(H) is fully self- 
consistent, and the value of the superconducting order 
parameter can be found from it for any temperature and 
geometry (i.e. with a proper account of the proximity ef- 
fects, where applicable). The self-consistency condition 
would easily follow from the variation of the action (|l3| ) 
with respect to A and finding the optimal configuration 
for the fields. However, it is convenient to impose the 
self-consistency requirement only at the very end of cal- 
culations. Any physical observable is then to be found 
by calculating an appropriate functional average with the 
functional (|13|)-(15). 

We proceed with illustrating how the model reproduces 
basic fundamental results for dirty superconductors, then 
demonstrate how to include consistently weak localiza- 
tion corrections in the vicinity of the superconducting 
transition in the presence of a magnetic field, and finally 
show how to take into account the self-consistency of the 
order parameter in the description of the proximity effect 
in the SNS geometry. 



IV. THE SIMPLEST APPROXIMATION 

We show that the basic results for dirty superconduc- 
tors can be reproduced in the simplest approximation: (i) 
we neglect all nonzero Matsubara harmonics of the pair- 
ing field, i.e. substitute A <5 ej _ e ' for A ee /; (ii) we neglect 
disorder-induced fluctuations near the saddle point, i.e. 
substitute the saddle-point value Q = A e . In this case, 
the matrix e+ A reduces to direct product over all integer 
n of (e„ + A ) ® (e n — A ) where 



A 



(16) 



Here Ao = |A| e zx is a two-component field which, natu- 
rally, plays the role of the order parameter (we omit the 
index in |A|). Now it is easy to find explicitly the eigen- 
values A and the diagonalising matrix Va in Eq. (nfl) . 



A e = \fei+ |A| 2 sgne„ , cos6> £ = e„/A e 



(17) 



V nA (r) 



cos — + 5 sgn e„ sin — 



where 5 = (Ao/|A|)(S ei _ e / is the 4x4 matrix which de- 
pends only on the phase x °f the field Ao and repeats 
the matrix structure of Ao, Eq. (|5|), and the full matrix 
Va is the direct product of all V n \ ■ 

On utilizing the assumption (ii) above, i.e. Q = A, and 
substituting the parameterisation ( |l7| ) into the Eq. (13), 



we arrive at the action S = Jd d r C with 



£ = 



SC = 



X T 
■nvD 



E 



(W e ) 2 



sin 



(18) 



Using the parameterization (|T^) one can easily sum over 
e to get 



6L = ^\c 1 (V\A\f + cJvx-^A) } cm, 



where the stiffness coefficients C\ 2 ar e given by 



LAT tanh 2T 



IAI 



1 

2T 



■ cosh 



2T 



(20) 



C 2 = 2|A|tanh^i. 



The functional (|19[)-(|2C|) coincides with that obtained in 
Rcf. [| Expanding coefficients C12 in A, one obtains the 
Ginzburg-Landau functional as that in Ref. |^. However, 
the simplest approximation used herejiand equivalent to 
those on which earlier considerationaJa were based) is 
not sufficient even in describing the vicinity of the super- 
conducting transition. In general, one must keep all the 
Matsubara components of the pairing fields. In the fol- 
lowing section, we will show how to do this in the vicinity 
of the transition in the weak disorder limit. 



V. GINZBURG-LANDAU FUNCTIONAL 

In the vicinity of the superconducting transition one 
can expand the action (|l3|) in the pairing field. A further 
simplification is possible by in the weak disorder limit, 
pp£ S> 1: one can integrate out the Q-field to obtain an 
effective action for the A-field only. In the quadratic in 
A approximation, the kernel of this action will give an 
effective matrix propagator of the pairing field, with due 
account for the disorder, which governs properties of a 
disordered superconducting sample near the transition. 

To integrate over the Q-field, one splits the action (|l? 
into S = Sq + Sa where 



5 = 



tivD 



Tr (d Q) 2 - —TreQ, 



(21) 



is the standard nonlinear a model functional as in the 
metallic phase. Then one makes a cumulant expansion, 
i.e. first expands e~( So+SA * > in powers of iSa, then per- 
forms the functional averaging with e _lS ° (denoted below 
by (. . .)q) and finally re-exponentiates the results. The 
expansion involves only the first and second order cumu- 
lants since the higher order cumulants generate terms of 
higher order in A. Then the only terms which contribute 
to the action quadratic in A are given by 



S off [A] = ^£ |dr|AJ 2 -^(Tr(A-e)Q 



^Tr[A A ,Q] 2 + ^(TrQ9oQA A f 



(22) 



4 



Expanding A and Aa to the lowest power in A and per- 
forming a standard functional averaging, as described in 
Appendix, one finds the action quadratic in A as follows: 



5 cff [A] = dr A^(r)(r| t u |r')A w (r') , 

with the operator IC^ given by 



(23) 



e(uj-e)<0 I 



i nf 2£ _, l (Q)c ' 

7Ti/ (2e-w) 2 



(24) 



Here nf^j(r, r') = ^r|n c,<i |r') are the cooperon and dif- 
fuson propagators, respectively, with 



C + \w\ 



(25) 



where the operator C = — D (V — 2ieA/c) defines the 
propagation of the cooperon modes; IT d is obtained from 
n c by putting A = 0. 

In the last term in Eq. @, Ilj^O) = II^|(r,r); this 
term may be obtained by expanding (in the weak disorder 
parameter) the cooperon propagator with the renormal- 
ized diffusion coefficient, 



C 



1 



1-^,(0) 



c. 



Therefore, this is just a weak localisation correction to 
the free cooperon propagator IIP ,(r, r'). 

The summation over Matsubara frequencies in Eq. ( [24] ) 
is easily performed to yield 



AttT 



(26) 



where T = T c0 (B = 0) is the transition temperature of 
the clean superconductor in the absence of a magnetic 
field. The weak localisation correction is proportional to 
the coefficient a u given by 



A-kT 



— T— 

lA/ ^ Dn 



Dq 2 



rn-Y ^ Dq 2 1 ' 



1 2 4ttT 



f 



4ttT 



1 2 4vrT 



For lo — the coefficient ao = a w= o(T') can be simplified 
in the two limits: 



no 



■nvL d 



E 



i 



Dq- 



V"(l/2) 



(27) 



8ir 2 vL d T ' 

where Lt = \/D/T is the thermal smearing length. 



The instability of the normal state (i.e. a transition 
into the superconducting state) occurs when the lowest 
eigenvalue of the operator K,^ becomes negative. The 
eigenf unctions of this operator coincide with the eigen- 
functions of the cooperon operator C. The lowest eigen- 
value of C is known to be Cq = DB/4>q, where 4>o is the 
flux quanta. This ground state cooperon eigenfunction 
corresponds to the lowest eigenvalue /Co of the operator 
ZC W . The condition /Co = implicitly defines the line 
T C (B) in the (T, i?)-plane where the transition occurs: 



In- 



Co 
47rT n 



(28) 



The term in the r.h.s. of Eq. ([28() describes a 1/g- 
correction to the main result. This weak localisation is 
linear in the magnetic field B and vanishes as B — > 
as expected (Anderson theorem). In a nonzero magnetic 
field the weak localisation correction to the B c is pos- 
itive which has a very simple explanation. The super- 
conductivity is destroyed by the magnetic field when the 
flux over the area with the linear size of the order of the 
coherence length becomes greater than the flux quanta. 
The weak localisation corrections diminish the diffusion 
coefficient which leads to the shrinkage of the coherence 
length. Therefore, one needs more strong field to ful- 
fill the condition of the coherence destruction. The same 
reasoning explains the growth of T c in the fixed magnetic 
field. 

Note finally that we have calculated the Q-averages in 
Eq. (HI) perturbatively, up to the first order in the weak- 
localization correction. It would be straightforward to in- 
clude main weak-localization corrections in all orders by 
calculating these averages via the renormalization group. 
This would lead to the renormalizing the diffusion coeffi- 
cient in the cooperon propagator (|25|) , thus changing the 
shape of the T C (B) curve. However, the value of T c (0) 
will again remain unaffected, since the superconducting 
instability is defined by the appearance of the zero mode 
in the operator K., Eq. (p6|). This zero mode is homo- 
geneous, and thus does not depend on the value of the 
diffusion coefficient in the cooperon propagator. 



VI. PROXIMITY EFFECT 

A recent supersymmetric versionS of the NLerM has 
been specifically formulated for studying the proximity 
effect in SNS junctions. Although this version is very 
convenient for a non-perturbative analysis, it has a nat- 
ural disadvantage of the supersymmetric approach: no 
interaction can be included beyond the mean-field ap- 
proximation. It means that the superconducting order 
parameter A should be treated as a background field 
rather than a dynamical one. More specifically, A was 
taken into accounted just by the boundary conditions 
(Andreev reflection) at the boundaries of a normal metal, 
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while having been considered as a given field in the su- 
perconducting region. This allows for changes in char- 
acteristics of the normal metal in the proximity of the 
superconductor, but not for the possibility of changes in 
the superconducting order parameter in the proximity of 
the normal metal. 

The_aetion in the normal region (N) has the standard 
forml§t£l while in the superconductor (S) we have NLctM 
of the form Eq. (|l3|) . The continuity of the Green function 
across the N/S boundary requires 



IN 



N/S 



Vl QsV A 



N/S 



(29) 



The N region by itself would favour Qn — A. The prox- 
imity leads to rotation of matrix Qn in the N region in 
order to match the structure imposed by the boundary 
condition (B9|) 



Q 



N 



Vl Q N Vi 



N 



(30) 



Xt\ N = x\ S ' cos6 £ \ 



N 



(34) 



Here |A| and x are modulus and phase of the order pa- 
rameter at the N/S interface. 

The classical trajectory corresponding-, to the action 
([52|) is nothing but the Usadel equation:E2l 



— A9 + - siir 2 1 
2 4 



2e , 

Vx A + ecos6> = 0, 



win- <y ( v x - ^A 



c 

= 0. 



(35) 



For quasi-l-D geometry in the absence of a magnetic 
field, the Usadel equation (|5|) can be written as the equa- 
tion for 9 



with the rotation matrix Vn of the same structure as 
Va in the S region so that at the boundary they match 
each other. Proceeding in the same manner as above we 
keep only the u = component of the pairing field and 
neglect the disorder induced fluctuations, i. e. we put 
Qn = Qs = A. Then for the N region we have 



VlAV N = cos C + sin £ 5, 



(cosx £ f^ r -l-sinxeff) 



(31) 



i(J 2 



where 6 t and \e are now independent variables. In a bulk 
superconductor, all these parameters were explicit func- 
tions of A and e, Eq. (|l7|). There is no such a constraint 
in the normal region. The (e, — e) sectors in the normal 
region are still coupled due to proximity effect but they 
may all be different. 

In this approximation the action corresponding to the 
iV region decouples into the sum of uncorrelated contri- 
butions: 



iSjv = 



(32) 



£e = 



D 



(V6 e y + sin" 9 e v x 



2e 



A 



ecost 



Now we find supercurrent j s by varying the action (|32| ) 
with respect to the vector potential A: 



j s = 2eiwD T (sin 2 # e (v X . 



^A 

c 



N 



(33) 



where (...)n stand for the functional averaging with 
the action Eq. (|32"|) , the functional integration being per- 
formed over the functions obeying the boundary condi- 
tions 



d 2 8 9 cos6> 9 

71 + Q ^ + 2 sm0 = O, 

ax z sm 6 



with the self-consistency condition on a e 



(36) 



da: 

XN=a e I — — 
sm v 



(37) 



Here \n = x+ - X 
two superconducting banks and L 



is the phase difference between 
^D/2e is the 



coherence length for two particles with the energy dif- 
ference e propagating in the normal metal. For a long 
normal bridge between the two superconducting banks, 
L ^> Lt = \/D/2nT, one may consider separately three 
regions: those close to the N/S boundaries (with the 
width of order Lt) and the bulk. Matching the solutions 
for all the regions, we find the following expression which 
well approximates the solution for the entire normal re- 
gion: 




FIG. 1. A spatial dependence of the phase \ E across the 
SNS contact for quasi- ID geometry. 
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6{x) = 8tan(^ /4) e - L/2i s/cos 2 ^ + sinh 2 — , 

V 2 L c (38) 

a c = 32tan 2 (# /4) sinxw Lj 1 exp [—L/L e ] , 

where 9o = 6 e \ N ^ s . In calculating the supercurrent 
through the normal bridge, one reduces the expression 
within the angular brackets in ( [33] ) to sin XNOt e . Then it 
is enough to keep only the leading term with eo = ttT 
because the contributions from all the other frequencies 
are exponentially suppressed as L e < Lt- Then we ob- 
tain the following expression typical for the Josephson 
junctions j s = j c syiixn, where j c is the critical current: 

j c = e2 7 nvD T tan 2 (6» T /4) L T X exp [-L/L T ] , (39) 

with 6t = &e - 

The supercurrent in the superconducting banks is 
found by varying the action ([ii]) valid in the S region 
with respect to the vector potential A: 



j s = eTTZADlAltanh^J- — , 
J II 2T L s 



(40) 



where is the length of the superconductor and xs the 
phase difference between its edges. 

It should be stressed that we have varied the action for 
the entire SNS structure, rather than only for the normal 
region subject to the boundary conditions at the super- 
conducting banks as in the supersypmetric variant of 
the NLctM for dirty superconductors£5l This means that 
the the phase difference across the normal region is not 
fixed but should be found self-consistently by finding the 
optimal configuration for the action for the entire SNS 
structure subject by the matching the fields at the N/S 
boundaries. This defines the actual phase difference xn, 
Eq. (|37]), across the normal bridge. Numerically, a sim- 
ilar procedure has been employed in Ref. |l]. It is easy 
to show that the matching condition can be expressed as 
the continuity of the supercurrents (as varying with re- 
spect to the phase difference is equivalent to varying with 
respect to the vector potential). Thus the supercurrent 
conservation defines the phase difference on the normal 
bridge: 



A , A Ls tit ■ 2 

-!— !-tanhJ— l -xs = —e~ L/LT sin ya, tan 2 — . 
64T 2T A L T A 4 



(41) 



so that if the width of superconductor banks L s is suf- 
ficiently large, the overall phase drop mainly happens 
across the banks. 

Finally, let us reiterate that the main result of the pa- 
per is a novel variant of the NLcrM given by Eqs. (jlg)- 
@. Here we have applied this formalism to a few rel- 
atively simple problems mainly to show that it works 
and has certain advantages over alternative variants of 
the NLcrM. This model has also been applied to a micro- 
scopic considerationE3 of the ctaaBtum phase slip problem 
in quasi-lD superconductors ,c3E3 and to a microscopic 



derivation of level statistics in nonstandard symmetry 
classes introduced in Ref. |25|. Let us also stressed that 
the method employed in the derivation of Eqs. (13)-jl5|) 
can be straightforwardly generalized both to including 
different types of interactions, and to considering the un- 
conventional pairing in dirty superconductors. 
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APPENDIX: 

To perform the functional averaging in Eq. (|2^), one 
should employ some parameterization of the-ficld Q in 
terms of unconstrained matrices, for exampleEilllj 

Q = (1 - W/2)A(1 + W/2)- 1 , 

where W = -W*, and WA+AW = 0. The Q-integration 
then reduces to the Gaussian one with the weight e~ s ° 
with So obtained from Eq. ( pl| ) by expanding Q to the 
second order in W. The Gaussian W integration is car- 
ried out with the help of the following contraction rules 



TrM W(r) P W(r') ) = x 



(Al) 



3/3/3 



(A2) 



ee' < 

a,/3 



(TrMVL(r)TrPVL(r')) = X 

ee' < 

a, (3 

where the upper indices a, (3 refer to the time-reversal 
sector and tr refers only to the matrix indices which are 
not i ndic ated explicitly. The matrix tt in the Eqs.(Al) 
and (A2) has the following structure in the time-reversal 
sector: 



7r ee /(r,r') = 



nf e _ e ,|(r,r') nf e _ e/| (ry) 



IT 



6-6' | 



,,(r',r) nji ,,(r,r') 



(A3) 



where the propagators are solutions to the standard 
cooperon and diffuson equations: 

[-DVl+oj] Ilj(r,r') = <S(r-r'), 
2e 



i — -A(r)V 

c 



(A4) 
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Note that in the absence of magnetic field these con- 
traction rules go over .to those previously derived for the 
orthogonal symmetry.Ea 

Next, one expands Q in Eq. ( ^20 up to the 4th power 
in W and uses the above contraction rules to obtain 



(Tr(A 



Tr(A — e)A . 



Tr [Aa, Q] ) = Tr [A a, A]" 

IQ 

nf 2e _ w ,(o) 



E 



e(w-£)<0 



(2e 



(V - —A) A, 



(A5) 



(A6) 



Taking into account that the second term in the brackets 
in the Eq. (|22| ) contributes to the higher order correc- 
tion only we arrive at the Ginzburg-Landau functional 
described in the text. 
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